LHCD модуль
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
type(Spectrum), | public | :: | full_spectrum | ||||
type(Spectrum), | public | :: | pos_spectr | ||||
type(Spectrum), | public | :: | neg_spectr | ||||
real(kind=wp), | public, | dimension(:), allocatable | :: | vvj | |||
real(kind=wp), | public, | dimension(:), allocatable | :: | vdfj | |||
integer, | public, | parameter | :: | kpt1 | = | 20 | |
integer, | public, | parameter | :: | kpt3 | = | 20 |
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!sav2008 pn=fn(r)
pn=fn1(r,fnr)
pn=fn2(r,fnr,fnrr) !sav2008
!!!!!!!read data !!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!sign of driven current in right coordinate system {dro,dteta,dfi}:
!!!curdir=+1.0 for current drive in positive direction "dfi"
!!!curdir=-1.0 for current drive in negative direction "dfi"
!!!spectrum Nz>0 is along dfi>0 and Nz<0 is along dfi<0
!!!it is also OK if Npar is used instead of Nz, but for Btor>0, that is along dfi>0
curdir=-dble(ispectr)
!!!!!!!!!!!!! begin iterations !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
find nevyazka
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(Spectrum) | :: | spectr | ||||
real(kind=8), | dimension(*) | :: | outpe | |||
real(kind=8) | :: | pe_out |
v12<v<=v2
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
double precision | :: | v1 | ||||
double precision | :: | v2 | ||||
double precision | :: | vmax | ||||
double precision | :: | cdel | ||||
integer | :: | ni1 | ||||
integer | :: | ni2 | ||||
integer | :: | ipt1 | ||||
integer | :: | kpt3 | ||||
double precision | :: | vrj(*) |
recalculate f' for a new mesh ! if(vrj(i).gt.vvj(i0)) exit
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | ispectr |
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real | :: | d | ||||
real | :: | u | ||||
integer, | intent(in) | :: | j | |||
integer, | intent(in) | :: | kmax | |||
real | :: | g |